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ABSTRACT 

The Elliptical Isolated X-ray (EllXr) Galaxy Survey is a volume-limited (< llOMpc) study of 
optically selected, isolated, L* elliptical galaxies, to provide an X-ray census of galaxy-scale (virial 
mass, Mvir^ objects, and identify candidates for detailed hydrostatic mass modelling. In 

this paper, we present a Chandra and XMM study of one such candidate, NGC 1521, and constrain its 
distribution of dark and baryonic matter. We find a morphologically relaxed hot gas halo, extending 
almost to R500, that is well described by hydrostatic models similar to the benchmark, baryonically 
closed. Milky Way-mass elliptical galaxy NGC 720. We obtain good constraints on the enclosed 
gravitating mass (M5oo=[3.8 ± 1.0] x IO^'^Mq, slightly higher than NGC 720), and baryon fraction 
(fb. 500=0. 13 ± 0.03). We confirm at 8.2-(t the presence of a dark matter (DM) halo consistent with 
ACDM. Assuming a Navarro-Frenk- White DM profile, our self-consistent, physical model enables 
meaningful constraints beyond R500, revealing that most of the baryons are in the hot gas. Within 
the virial radius, fb is consistent with the Cosmic mean, suggesting that the predicted massive, quasi- 
hydrostatic gas halos may be more common than previously thought. We confirm that the DM and 
stars conspire to produce an approximately powerlaw total mass profile {ptot cc r~") that follows 
the recently discovered scaling relation between a and optical effective radius. Our conclusions are 
insensitive to modest, observationally motivated, deviations from hydrostatic equilibrium. Finally, 
after correcting for the enclosed gas fraction, the entropy profile is close to the self-similar prediction 
of gravitational structure formation simulations, as observed in massive galaxy clusters. 
Subject headings: dark matter — Xrays: galaxies — galaxies: elliptical and lenticular, cD — galaxies: 
ISM — galaxies: formation — galaxies: individual (NGC1521) — galaxies: funda- 
mental parameters 



1. INTRODUCTION 



The distribution of mass in galaxies, both in the form 
of baryons and dark matter (DM) , is a crucial yardstick 
for elucidating galaxy formation and evolution. Dissipa- 
tionless DM simulations in our current (ACDM) cosmo- 
logical paradigr n predict DM halos with a characteristic 



extend to a lower mass regime the hydrostatic X-ray tech- 
niques widely employed t o study the dark and baryo nic 
matter in galaxy clusters (Buote & Humphrey 2012a, for 
a review) . Provided the hot gas halo is sufficiently bright 
and morphologically re laxed, hvdrostatic techniques are 



^00l| ; Maccib et al. ^008 ). In the centres of these ha- 
los, baryons condense into stars, but the complex inter- 



expected to be useful ( Buote fc Tsai _ 1995 ;_ Rasia ct al 



density profile (iNavarro et all & B. the average 2006|: [Nagai et alj 2007], Piftarm^ Valdarmm| gOOg; 
shape of which vari es slowlv VifFlniil Bullock et al.| pH^te Hump;E^Jp012a|) . This is confirmed observa- 



Daryons condense into stars, but the complex inter- 
play of gas cooling and heating involved (including ac- 
cretion shocks, feedback from supernovae, stellar winds, 
and 



ictivc galactic nuclei) has been th e aubicct of vig ' 



orouil debate (e.g. [Whito fc Rccs||1978|; [White fc Frenk 



1991 



2006) 



Kcrcs et al. 20O5 ' Hopkins et "aT ' 2t)(7^ Croton ct al 
. Current models tend to predict that large {^L^,} 
galaxies should possess massive, but diffuse (a nd, possi- 
bly, hard to detect) coronae of hot baryons (e.g. Fukugita 
|fc Peebles||200^ ; [Grain et al.||2010| ). 

Giant elliptical galaxies provide a natural laboratory 
for e xploring those predictions. The X ray omitting 



hot paw Vialns; aynnnr l many par1y-typp galavipH (p g 

'Sullivan et al. 2001) provide a unique opportunity to 
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Buote &: HumphreyT 2012a ). This is confirmed observa- 
tionally by detailed comparisons between masses inferred 
from X-ray methods and independent stellar dynamical 
measurements, or the predictions of stellar population 



synthesis models, whic h suggest a ty 
at lea s t ~20 - -30% ( e.g. IChurazov et al 



jical accuracy of 



et al.[ b008[ . [2009a[: IShen fc Gebha^pOTo| ; |Das et al 



20To| ; poTl| [Humphrey et al.[[2012b| )" 



2008 : [Humohrev 



!STthough the properties of the largest galaxies may be 
inter twined with that of a surrounding group or cluster 
(e.g. [Helsdon et"aL[ p001| [Mathews et al.[[2006| ), there 
is increasing evidence that ^L* early- type galaxies in 
~Milky Way-mass halos can be found with hot gas de- 
tectable out at least to ~R25oo- In particular, Humphrey 



et al. ( 2006 ) presented a hydrostatic analysis of three iso 
lated galaxies, based on data from the Chandra X-ray 
Observatory, inferring virial masses ( Mvir) < lO-'^'^Mr:^ 



Using, deep Chandra and Suzaku data, [Humphrey et al" 



(2011) refined the analysis for one of these systems, 
NGC 720, confirming Mvi,^( 3.1 ±0.4) x 10^-^M,t._, close to 



the mass of the Milky Way ( [Klypin et al.|[2002D . Other 
X-ray bright, fairly isolate d galaxies with similar proper 



ties have been discussed by|0'Sullivan fc Ponman (2004), 
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p 'Sullivan et alj ( p007l ) and [Memola et al| ( 1201 

'i'he extent to whicn typical giant elliptical galaxies 
possess DM halos i n accord with ACDM remains unclear. 
p3uote et al. ( 2007 ) assembled from the literature hydro- 
static measurements of the virial mass and concentration 
(Cvir=Rvir/?'s, whcrc Rvir IS the virial radius and is 
the characteristic scale of the DM halo density profile) 
for a sample of galaxies, groups and c lusters, including 
t he th ree isolated systems studied by Humphrey et al.| 
(2006). They found, for the first time, an mverse cor- 
relation between Cvir and Mvir, with a slope close to the 
theoretical value, confirming a fundamental prediction of 
ACDM models. Still, the three isolated galaxies, consti- 
tuting the low mass (< lO^^M©) end of the relation, were 
marginally ('^2.6-cr) more concentrated than predicted 
by recent theoretical models in the most favourable cos- 
mo logy (the "WMAPl " model for relaxed halos reported 
by Maccio et al, _ 200g ]^ Given the small number of ob- 



ally lack the resolu tion to resolv e the mass profiles in 
detail. Indeed, the Auger et al. data may be slightly 



offset from, and exhibit more intrinsic scatter than, the 
X-ray relation, although more X-ray data are needed. 

No census of the mass within early-type galaxies is 
complete without accounting for the baryons in the hot, 
diffuse gas. In the local Universe, the measured stellar 
and cold gas content of galaxies li es significantly 



below the Cosmologi cal baryon fraction ( Fukugita et al, 
19981 ; iMcGaugh et al J ^oTo|) , ' and m tension with standard 



models of galaxy formation (e.g. Benson et al, 2005| ) 



Resolutions to this problem generally involve the "miss- 
ing barvons" eith er resi ding in a massive, h ot halo (e.g 
Mailer k Bullock| ^004| ; [Fukugita fc Peebled |2006|) or be- 



jects involved, however, it is unclear whether this repre 



sents 



ing e j ected completelv from the sy stem (e.g. |Dekcl fc Silk 
1986 ; Oppenheimer &: Dave 2006 ). To date most of the 
effort to locate these hot halos has focused on disk g alax- 



actual tension with theory, or is a consequence of 



selection ettects or small number statistics. 
To 



date, these represent arguably the best constraints 



ies. where they have not yet been robustly de t ected ( Ben 
son et al. 2000| ; Anderson & Bregman 2010; Rasmussen 
et al.| 2009). Although X-ray absorption line stu dies have 
identified hot gas around t h e Milky Way (e.g. [ Nicastro 



on t. le DM halo concentration tor individual giant el 



liptical galaxies that reside in galaxy-scale {< 10^'^Mq j 
halos. Lensing studies alone cannot presently resolve 
the mass profiles of individual galaxies, although stacked 
weak (and strong) lensing analysis indicates that at least 
some early-type galaxies can be foimd in lO^^M^ 



et all ^002|- [Fang et al.|[2002|: [Rasmussen et al 

et al.l^nnS ; Bregman i^FTToyd-Davies [2007 ~ ljuote' et al 
20091; Fang et al.[[2010 ) , and there have been reports of ex- 



2003[: [Fang 



DM jialoes with concentrations broadly consistent with 



ACDM ( [Mandelbaum et al.[ B006[; pavazzi et al.[ E007[). 
Orbit- and particle- based stellar dynamical methods are 
beginning to emerg e that incorporatj s DM halos (e.g 



tended X-ray emi ssion unassociated with star formation 
in a disk galaxy (Anderson & Bregman 2011), whether 
these constitute the predicted major reservoirs of ba ryons 
depends on, generally uncertain, extrapolation (e.g. Fang 
et alj 200£j Rasmussen et al. [2009 ; Anderson &: Bregman 



[Thomas et al. 2007 ; de Lorenzi et al. 2009), but there 



have, so far, been few published constraints on c 

If the DM halos o f early-type galaxie s are well de- 
scribed by the NFW (Navarro et al. 1997) profile, within 
~the optical effective ra dius (R^), tne ba ryonic compo- 
nent must be dominant (Buote et al. 2012| ). It is increas- 
ingly being recognized that the baryons and DM conspire 
to produce a total mass density profile that can be well- 
approximated by a powerlaw (pt nt oc r~") over a wide ra- 



dial range (e.g. [ Fukazawa et al. 2006; Gavazzi et al. 2007 
Humphrey fc Buo^ _ 2010 , and refe rences therein: Uhu , 
razo ^ et al.][201(j| ; [DuSyct al.[[20To[) Although typically 



20T0[ , p0Tl[ ) 

Recently Humphrey et al.[ ( ^011 ) studied the hot gas 
around the isolated, ^Milky Way-mass (Mvir =[3.1 ± 
0.4] X 10^2) elliptical galaxy NGC 720, detecting the 
baryons as far as ^ i?2500- Elliptical galaxies have a 
distinct advantage over disk galaxies for detecting the 
putative hot halos, since few of the baryons are bound 
up in cold gas, so the halo should be denser and more 
luminous than in a comparable spiral galaxy. In the 
case of NGC 720, the baryon fraction within R25oo0 was 
tightly constrained to fb,2500= 0.10 ± 0.01, rising to 
fb= 0.16 ± 0.04 by Ryj,,, consistent wi th the Cosmolog- 



ical value (0.17: Dunkley et al. 2009). Unlike the ad 



2 i s repo rted (Koopmans et al. 2009), Humphrey & 
Buot^ ( [2010 ) found that a derived irom Chandra obser- 
vations of a sample of 10 galaxies, groups and clusters 
was tightly anti-correlated with Re. This behaviour can 
be understood by the combination of the stellar (Sersic) 
and DM (NFW) profiles required to maintain an approx- 
imately powerlaw total mass distribution, and implies 
that the DM fraction within Re varies systematically in 
such a way as to reproduce, without fine tuning, the tilt 



hoc extrapolations often employed (e.g. the isothermal 
/3- model) , this self-consistent evaluation of the model at 
large radii required only that the gas is approximately 
hydrostatic (as exp ected around an isolated system, e.g. 



Crain et al 
to N_b'W, w 



2010|), and the DM mass profile is close 
lie being relatively insensitive to the ther- 



modynamics of the gas outside the region where it was 
clearly detected. NGC 720 therefore constitutes the most 
promising detection of a baryonically closed^ Milky Way- 
mass halo, which indicates that feedback need not denude 



of the fundamental plane (FP). If this trend is the ex- 



such a galaxy of a large fraction of its baryons (e.g. Kauf 



planation of the tilt of the FP, it should exhibit very 
little intrinsic scatter, but more X-ray measurements are 
needed to investigate this further. The anti-corr elation 



betwee n a a nd Ro has recently been confirmed by Augei 
et al[| (2010), who used joint strong lensing and stellar 



kinematics measurements. While the optical results pro- 
vide important verification of the X-ray work, they are 
limited by being confined mostly to within ~Ro, where 
the DM halo is sub-dominant, and because they gener- 



mann et al. 2009). Whether or not NGC 720 represents 
an unusual case remains to be established, but may have 
clear implications for the location of the missing baryons 
in the local universe. 

Although X-ray observations of galaxy-mass 
(<101^Mq) halos provide a unique and powerful 
insight into galaxy formation, the small number of 

^ We define as the three dimensional radius within which 
the mean mass density of the system is A times the critical density 
of the Universe. 
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reliable measurements of Cv 



baryon fraction (fb), 



and a at this mass scale limits the conclusions which can 
be drawn. In addition to the problem of small number 
statistics, the current objects were chosen for study het- 
erogene ously, potentially introduci ng unknown selection 
effects ( Buote fc Ifumphrcy 2012a , for a review). To ad- 
dress these concerns, we initiated the Elliptical Isolated 
X-ray (EllXr) Galaxy Survey, an X-ray survey of an 
optically selected, volume limited (< 110 Mpc) sample 
of very isolated, '^L* early-type galaxies. The isolation 
condition eliminates most group-scale halos and ensures 
an accurate census of the X-ray properties in systems 
resembling the optical properties of NGC 720. Isolation 
also minimizes the likelihood that the low-density hot 
gas halo expected around a ~Milky Way-mass galaxy 
is stripped in a dense, cluster environment. A full 
description of the sample, and initial results are given in 



Buote et all (pOll). 

^ni> 



INot ail ot the EllXr galaxies have luminous X-ray halos 
within ~Rc, as might be expected given the large range 
of measu red Lx at fixed optical l uminosity in early-type 



galaxies ( Canizarcs ct al 
Elhs fc q^iiivanl [jUU6| ).' 



1987| ; pS 
I'his largi 



Sullivan et al.| |2001 

ge scatter may indi- 
cate that both the virial mass of the halo and the feed- 
back history of the ga laxy play a role in determining Lx 



(Mathews et al. 2006). Since the gas emissivity depends 



on the square o f its density, and gas density profiles are 



1 



not self-similar (Humphrey et al 
" 007|; |Sun et"aL| 200H|), it is difficu 



2006: Gastaldello et al 



t to map galaxies from 
the optical versus X-ray luminosity plane onto fb directly. 
This means that Lx, which is typically measured within 
only a small fraction of the virial rad ius (Rvir), is a poo r 



tracer of the overall gas mass (e.g. Grain et al. 2010). 



While even low-Lx early-type galaxies could possess mas 
sive hot coronae in which most of the gas has been pushed 
out to large scales, detecting them will pose si milar obser- 
vatio nal ch allenges as for spiral galaxies (e.g. Rasmussen 
et alj |200S ). Instead, those objects with X-ray luminous 
halos provide arguably the best opportunity to measure 
the baryons out to '^tens of kpc, and thus provide direct 
constraints on the baryon fraction, at least for a subset of 
early-type galaxies. It is therefore of interest to define an 
initial X-ray luminous sub-sample of the EllXr galaxies 
for further study with detailed hydrostatic methods. 

In this paper, we present a detailed Chandra and XMM 
study of one such galaxy. NGC 1521 was identified as 
one of several X-ray luminous analogues to NGC 720 in 
a shallow, pointed XMM observation taken as part of 
the EllXr Galaxy Survey, and was targetted for deeper 



follow-up. 
detai l in § 



We discuss the properties of the galaxy in 
^, bufoic dcsciibing the X-ia^i dAtA-ieductiun 
and analysis (§ the mass modelling method (§ U), 
the likely sources of systematic uncertainty (§ |^) and 
reaching our conclusions in § ^ 

We adopted a distance of 61.2 M pc to NGC 1521, cor - 
responding to the redshift 0.01415 ( Ogando et al. 2008 ), 
if we assume a flat cosmology with Hq = 70km s^^ and 
r^A = 0.7. In § m we show that small errors in our dis- 
tance estimate will not affect our conclusions. At that 
distance, 1" corresponds to 290 pc. We adopted i?io2 as 
the virial radius fR,Hr). based on the approximation of 
Norman] (p98|) for the redshift of NGC 1521. 



Bryan 



Unless otherwise stated, all error-bars represent l-cr con- 
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Fig. 1. — niatHKiitinn nf gnlgy^cg around NGC 1521 (shown as 
a large star) in the \HyperLEDJ^ database. Galaxies with known 
redshift are shown as large triangles (less than 2 apparent mag- 
nitudes fainter in the B-band than NGC 1521), large circles (2- 
4 magnitudes fainter), or small circles (more than 4 magnitudes 
fainter). Galaxies with measured recessional velocities differing 
by more than 1125km s~^ from NGC 1521 are omitted from the 
upper panel; none of the galaxies in the upper panel are within 
2 magnitudes of NGC 1521. The corresponding redshift range is 
shown as dashed blue lines in the lower plot. Those galaxies in the 
field of view but without redshift information are shown as large 
squares (2-4 magnitudes fainter than NGC 1521 in the B-bancLoL- 
— LL-unaiiailable, in the I-band, Ks or "Opt" band, as given in \Hy-\ 
perLEDA , respectively) , or small squares (more than 4 magnitudes 



tamter j. The red squares indicate the three galaxies catalogued, 
airing mi+h ^(^f^ 1 -'''^l as belonging to the small "group" S138 by 



Ftamella et al.| ( |2002| ). The dotted ellipse indicates the region en- 
closing projected K500, while the dashed circle corresponds to Rvir- 
These regions appear elliptical due to the aspect ratio of the figure. 
NGC 1521 is clearly isolated from other bright galaxies. 

fidence limits (which, for our Bayesian analysis, implies 
the marginalized region of parameter space within which 
the integrated probabflity is 68%). 



2.1. 



2. NGC 1521 
Target selection 



The EllXr galaxies were identified in the HyperLEDA 
database to be optically isolated, ~L* early-type galaxies 
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within 110 Mpc. For a full description of th e sample se- 



l ection and properties, we refer the reader to Buote et al 



(2012). Briefly, we selected early-type galaxies with ab- 
solute B-band magnitude between -21.4 and -19.8, and 
required there to be no other galaxy within a projected 
distance of 750 kpc that is less than 2 apparent magni- 
tudes fainter than the target. Where available, we used 
recessional velocity information to eliminate foreground 
or background interlopers (for which we assumed the line 
of sight velocities differed from the target by more than 
1125km s-i). 

In Fig ^ w e show the distribution of galaxies in the 



HyperLEDA database around NGC1521. Excluding an 
obvious, bright foreground object (NGC 1518), none of 
the galaxies are less than 2 magnitudes fainter than 
NGC 1521, confirming its isolation from bright compan- 
ions. We note t hat, n evertheless, the group catalogue of 
Ramella et al. (2002), derived using a friends-of- friends 
algorithm, mcluded NGC 1521 as a member of the low- 
mass (Mvir^i; IO^Mq) S138 "group". The other three 
members of the putative group are marked in Fig n|, but 
all are more than 2 magnitudes fainter than NGC 1521 
(so that it does not violate our isolation criterion). Fur- 
thermore, while these galaxies could be gravitationally 
bound to the NGC 1521 system (as the brightest mem- 
ber), they lie outside its virial radius (§ H), and are 
not symmetrically distributed about NGC 1521, implying 
that SI 38 does not represent a virialized system (whereas 
NGC 1521 itself does). 

2.2. Optical properties 

In the optical, NGC 152 1 exhibits complex st ructure. 
B-band isophotal analysis ( |Capaccioli et al. 198g ) reveals 
an ellipticity gradient and a strong position angle twist 
with in the central ^10 kpc. indicating that the galaxy 



is no |t axisymmetric, and may be triaxial. Although our 
isolation criterion ensures that NGC 1521 is not currently 
undergoing a major galaxy interaction, deep B, V, R and 
I images ta ken as part of th e Carnegie- Irvin e Galaxy Sur- 
vey (CCS: |Ho et al.| |201l| ; |Li et al.| ^01lD p| reveal faint 
shell-like surface brightness discontinuities outside the 
central ~40 kpc, which may be the relic of a past mi- 
nor merger event. Still, we note that such features may 
persi st in the stars for several Gyr following the event 
(e.g. Bchweizer & Seitzer |1992| ), while the sound cross- 
ing time ot a ^^40 kpc system, assuming an ambient gas 
temperature ^^0.5 keV/k (similar to what is measured 
), is only ^ 10^ yr. This means that 
relax to an approximately hydrostatic 
state following a significant disturbance far faster than 
the stellar distribution, and so it is unsurprising that the 
X-ray morphology appears very relaxed (§ ^). 



in NGC 1521; § 
the hot ISM wil 



Accurate modellmg ot the mass distribution requnes 
a refiable depiuJecLiuii of Lhe sLellai lighL piufile. — 
obtain this, we explored a triaxial model for the un 
derl y ing s tellar light distribution. Following [Cappel 
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Fig. 2. — Top panel: I-band isophotal major-axis position angle 
(measured counterclockwise from the north) of NGC 1521 (Hata 

pn;^tg^ rloriynH frnm flip r|arncgip-Tv^Hnp Galaxy SurvCy data 



ct al. 2011; Li et al. 2011). Overlaid is our best-fitting triaxial 
model (solid red Imc). Uenter panel: I-band ellipticity profile and 
best-fitting model. Lower panel: I-band major axis surface bright- 
ness profile of NGC 1521, shown along with our best-fitting model. 
The fit is good; the mean absolute residual is only 0.01 mag, which 



is sufficient for our p 
published profile of 



ivpricna 



Note-that fij differs slightly from the 
(2011) since we allow the isophotal 



Li et al 

position angle and ellipticity to vary with radius during its com 
putation. 

y and z are Cartesian coordinates. To simplify the de- 
projection, we assumed that the long, short and interme- 
diate axes of each ellipsoid were coaligned. For any given 
orientation (which is completely described by three posi- 
tion angles, 0, (j) and ip), we projected th e ellipso i d ont o 
Cartesian sky coordinates, using Eqn 9 of Binney ( |1985| ). 
We used dedicated software to fit this model to the cen- 
tral 4.3x4.3' portion of the calibrated , fiat-fielded I-ba nd 
image prod uced by the CCS team (Ho ct al. 2011; O 
et al. 2011), while masking unrelated point sources ana 
other galaxies. To account for seeing, we convolve the 
model image by a 1.2" FWHM Gaussian, estimated from 
the broadening of stellar images near NGC 15 21. We 



lari| dgOOg), we approximated the stellar density profile adopt ed the Galactic extinction correction from [Schlegel 



as a iseries oi concentric, triaxial ellipsoids with Gaus 
sian ladial piofiles (/j, ix exp {~0.b{uv /(jY)), wlieie 
is the stellar density, ay is the ellipsoidal coordinate 
ay = (x)^ + (y/p)^ + {z/qY (here I > p > q), and x. 



http: / /cgs. obs.carnegiescience.edu/CGS/Home.htmj 



et al. (199S). We allowed the position angles of the sys 
tem, and the luminosity, cr, p and q of each component to 
vary during the fit. We obtained a satisfactory fit when 
we employed 6 Gaussian components, and not only accu- 
rately recovered the surface brightness profile (as shown 
in Fig H), but also the isophotal position angle twist to 
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~ 0.5° and the ellipticity to ^ 2%, on average. 

3. X-RAY DATA REDUCTION AND ANALYSIS 
3.1. XMM 

The region of sky containing NGC 1521 was imaged 
by XMM on 4 separate occasions, as part of the EllXr 
Galaxy Survey. To simplify the analysis, we consider here 
only the results from ObsID 0552510101, the deep, 104 ks 
(excluding periods of background flaring) exposure be- 
ginning on 2009 Feb 7. In our default analysis, we focus 
only on the EPIC MOS instruments, as we have found 
the calibration of the PN to be more uncertain in the 
low-temperature (kT^O.5 keV) regime, and the source 
is too faint for the RGS data to be useful. Still, as a 
systematic error chec k, we explore the results obtained 
with the PN in § The data were reduced and anal- 
ysed as described m Gastaldello et al. ( [2007 ), using the 
XMMSAS 10.0 software suite. Briefly, calibrated events 
files were generated with the emchain task. We filtered 
the list to remove known bright pixels and hot columns, 
and periods of high background (flaring) were identified 
visually in the 10-12 keV lightcurve, and excised. Data 
from the unexposed portions of CCD MOSl-4 exhibited 
an unusually high count-rate and soft spectru m, consis- 
tent with the "a nomalous states" identified by Kuntz & 
^nowdeii (2008). We therefore excluded all data from 



this CCD in our subsequent analysis. Point sources were 
identified by operating on the 0.5-10 keV band images 
with the ewavelet task. The source lists from the three 
EPIC instruments were merged and checked visually. In 
subsequent analysis, we excluded data from within a 30" 
radius aperture (corresponding to approximately 90% en- 
circled energy) centred on each confirmed source. The 
image of the central MOSl CCD is shown in Fig |[ re- 
vealing both the extended emission from NGC 1521, and 
a number of point sources that were excluded. 

Spectra were accumulated in 8 concentric, contiguous 
annuli, placed at the X-ray centroid, reaching ~190 kpc. 
To mitigate mixing between annuli due to the point 
spread function, the minimum half-width of the an- 
nuli (corresponding to the radius of the central, circu- 
lar region) was set to 30", corresponding approximately 
to ~90% encircled energy. Typical spectra are shown 
in Fig 0. Spectral redistribution matrix files (RMFs) 
and ancillary response files (ARFs) were generated with 
the SAS rmfgen and arfgen tasks, the latter using an 
exposure-corrected detector map to perform fiux weight- 
ing. 

The MOSl and M0S2 spectra were fitted simul- 
taneously in Xspec vers. 12. 5. In, to obtain the pro- 
jected abundance, tem perature and density profiles (see 
Humphrey et al. 2011). The data were fitted using the 
C-statistic, which is less subject to bias in all count 
rate regimes than the pop ular implementations of for 
Poisson-distributed data ( [Humphrey et al" ^009b| ), and 
the fits were restricted to the 0.5-5.0 keV band (to avoid 
instrumental lines at higher energy). To aid convergence, 
we rebinned the spectra to ensure at least 20 photons per 
bin. We modelled the hot gas as an APEC model, and in- 
cluded a 7.3 keV bremsstrahlung component to account 
for undetected LMXBs, which was only sig nificant withi n 

Rci(4 7 kpr, from the 2MASS database- [jmrrti] [2000D 



tween all annuli. Where they coul d not be constrained , 
they were fixed at the Solar ratio ( Asplund et al, 2004 ) . 
To acco unt for the backg r ound , we adopted the ap- 



proach in Humphrey et al. (2011). Specifically, to ac 



count for the instrumental and particle background, we 
included a broken powerlaw model (not folded through 
the ARE) and two Gaussians (at 1.5 keV and 1.7 keV). 
The normalization of the instrumental components was 
allowed to vary with radius, and (to improve constraints) 
the shape of the broken powerlaw component was tied 
between the inner three annul i. and between the fourth 
and fifth annul i. We n ote that |Kuntz &: Snowden ( 200^ ); 
Bnowden et al. (2008) provide an alternative strategy for 
backgroun d su btraction (ESAS|), which we explore in 
detail in 5 5.4. We choose not to use ESAS for our de- 



fault analysis since there are inherent uncertainties in the 
procedure to map from the out of field of view count rates 
onto the instrumental background, which could cause 
problems in the highly background-dominated regime at 
the outskirts of NGC 1521. To account for the cosmic 
X-ray background, we included an (absorbed) powerlaw 



model with r=1.41 ( pe Luca fc Molendi] p00"4| ) . 

Given its Galactic coordinates (1=216", b=— 45°), 
NGC 1521 is located at the extreme edge of the larger 
of two adjacent excesses in the soft X-ray background, 
known as the "Eridanus X-ray enhancer nent" . that may 
be due to an old supernov a remnant ( Naranan et al 



1976; Bnowden et al. 1995). Therefore, it is impor 



taut to take care over the characterization of the soft 
X-rav background. Using pointed XMM observations, 
Henley et al ( 2010 ) were able to characterize the Galac- 
tic background at high latitudes, including data in the 
vicinity of NGC 1521, with a model comprising two 
plasma components (one unabsorbed and one absorbed) , 
with kT^O.l keV and ^ 0.2 keV, respectively (see also 
Kuntz & Snowden 200(]| ). For the two pointings within 
-8° of NGC 1521 (tEeir observations 22 and 23), kT 
of the hotter plasma component was constrained to 



0.197^ 



keV and 0.189 



+0.009 



-0.018 

are consistent within errors 



Q Qj^4 keV, respectively, which 



The ROSAT 0.1-1.0 keV 
count-rate in a 0.6-1.0° region around N GC 152l[| isac- 



tually within ~15% of that coincident with Henley et al 
observation 23. 

Since NGC 1521 is closer to the peak of the Eridanus 
enhancement than these fields, the temperature of the 
hotter component could, plausibly, be different. To ex- 
plore this, we analysed a deep (~100 ks) archival Suzaku 
observation of the "Eridanus Hole" (ObsID 502076010), 
a blank-sky region also in the outskirts of the Eridanus 
enhancement (1=213°, b=— 39°), but where the back- 
ground is ~30% higher than for NGC 15 2 1. We reduced 
the data as described in Humphrey et al. ( 2011 ) , and ex- 
tracted a spectrum from the whole field, excluding data 
in 2' regions around bright point sources (we identified 
these sources by eye in the XMM EPIC MOSl image 
of the same field (ObsID 0203900101); we used Suzaku, 
rather than XMM for characterizing the Galactic emis- 
sion due to its lower, stabler instrumental background 



^ http:/ /hcasarc. gsfc.nasa.gov/docs/xmm/xmmhp_xmmesas. html 



' JJctpaainrrl hy qnrryin 
maps ( Snowden ct al 

y-^ay harkgrniinri Tnn 



the. KU S AT diHusc X-ray background 

1997 ) . wHh thp nn-lir,o HFABARP 



Hot |as abundance ratios with respect to Fe were ticdTx;- bin/Tools/xraybg/xraybg.p: 



http; / /hcasarc. gsfc.nasa.gov/cgi- 
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Fig. 3. — XMM MOSl (central CCD; top left) and Chandra (top right) images of NGC 1521. Point sources found in the XMM image 
have been marked with white circles, and are excluded from subsequent analysis. Although it is not clear given the dynamic range of the 
JfMM image, the emission from NGC 1521 extends to the edge of the central CCD, as we show in §J3. The Chandra image has been cleaned 
of point sources and mildly smoothed, and arbitrarily spaced isophotes overlaid to guide the eye. The smoothing scale varied from ~1" at 
the smallest scales, to ~0.9' (15 kpc) at the outer part of the image (see text). To explore the slight distortion of the innermost isophote, 
we show (bottom left), the azimuthal variation of the surface brightness in the (unsmoothed) Chandra image, averaged radially between 
1-5" (upper, black data-points) and between 5-15" (lower, red data-points). The position angle is measured anticlockwise from due east. 
We overlay the best elliptical /3-model fit (allowing a position angle twist between each region). We note that the slight excess between 
240-300" in the 1-5" data is not statistically significant. At the bottom right, we show a "residual significance" image (see text) of the 
centre of the system, indicating deviations from a smooth fit to the Chandra X-ray isophotes. There is no obvious large-scale feature in 
this map, indicating the system is largely relaxed. 

and good spectral resolution). As this is effectively a 
blank-sky field, after subtracting off the standard in- 
strumental background (which was estimate d by stan- 
dard Heas oft tools; for more details, refer to Humphrey 
et al. 20f 1), the remaining spectrum represents the X-ray 
We were able to fit this in the 0.5-5.0 keV 



used the whole Galactic column density (Kalberla et al 
2005) to absorb the hotter component 



background 

band with a model comprising a F = 1.41 powerlaw, 
an unabsorbed .07 keV APEC plasm a model (with So- 
lar abundances: Asplund ct al. 2004), and an absorbed 
APEC model with kl'=0.22 ±Um keV, i.e. approxi- 



r natelv the same model as was used by Henley et al 
( 2010 ). Although the inferred distance of the gas con- 
trib uting to the Eridanu s enhancement is close (~400 pc; 
e.g. ^nowden et~aL||l995| ), we found an adequate fit using 
the whole Galactic column for the absorbed APEC com- 
ponent, which may reflect the high galactic latitude of the 
feature. Freeing Nh did not give rise to a significantly 
better fit. We therefore adopted kT of 0.07 and 0.20 keV, 
respectively for the two APEC components in our fit, and 



We found this background model was able to fit ad- 
equately the background spectra generated for regions 
corresponding to each annulus used in our XMM analy- 
sis from the standard "template" events files. We found 
that our full model was able to fit the spectra in all annuli 
reasonably well, as shown in Fig ^. We show the inferred 
Fe abundance profile and the global abundance ratios in 
Fig ^, and the measured temperature and density pro- 
files in Fig H We note that we were able to obtain good 
constraints on the gas temperature and density out to 
190 kpc, which is well beyond R2500J a-nd is approaching 
R500 (R2500=121±8kpc; R5oo=240±22 kpc; §1), which 
is beyond even what we achieved in our NGC 720 analy- 
sis. We expect the results not to be strongly sensitive to 
the treatment of the background if the ratio of the source 
to the background rate >0.5 in the 0.65-0.9 keV band, 
corresponding to the Fe L-shell region, which is the cru- 
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Fig. 4. — Representative Chandra and XMM MOSl spectra for NGC1521, shown without background subtraction. In addition to the 
data, we show the best-fitting model, folded through the instrumental response (solid black line), along with the decomposition of this 
model into its various components. We show the hot gas contribution (sold red line), the composite emission from X-ray binaries (dash-dot 
magenta line), the instrumental background (dotted orange) and the cosmic X-ray background (dashed purple line). The background is 
dominated by the instrumental component, but emission from the ~0.5 keV gas is detectable above the background below ~1 keV in all 
the spectra. 



cial temperature diagnostic for a ^^0.5 keV plasma. This 
is true for the inner 5 XMM annuh, but for the outer 
3 annuh, the source/ background ratio fahs to ~0.22, 
0.08 and 0.07, respectively. We explore the impact of 
diffe rent background subtraction subtraction methods in 



detail in § 5.4 , finding, unsurprisingly that the outermost 
annuli are most sensitive to this choice, as shown in -b'ig|^. 
Nevertheless, we do not reach qualitatively different con- 
clusions regarding the global parameters of the system 

(§|53)- 



The region of sky containing NGC 1521 was imaged 
by the ^ C/S instrument aboard Chandra (ObsID 10539; 
beginning on 2009 Jul 4) for a total of 49 ks good time 
(with periods of background flaring removed) . The data 
were r educed and analysed as described in Humphrey 



et"al| ( ^012aD , using the CI AO 4.3 software suite and 
the corresponding Chandra calibration database ( Caldb) 
vers. 4.4.2. Briefly, the data were reprocessed from the 
"level 1" events files, following the standard data reduc- 
tion threads^. Periods of high background were identi- 



3.2. Chandra 



http: / /cxc. harvard.edu/ciao/threads/index. html 




Ne Mg Si 
Element 

Fig. 5. — Upper panel: Projected Fe abundance profile, mea- 
sured with Chandra (triangles) and XMM (stars). Note the overall 
good agreement between Chandra and XMM. Lower panel: Aver- 
age abundance ratios with respect to Fe for the Chandra (triangles) 
and XMM (stars) data. The solid line is the best-fit model where 

the onrirlirr)(^nt rnmf^Q frgrri SNTo and SNII, aSSUmiug the W7 SNTg 

( 1997b) and the SNII yields from Nomotc 
enrichment fraction is 0.77 ± O.Ua. Fot 



yields, from D>[omoto ct al^ , 
et al.||( |l997a| ). 'I'he SNIa 



reference, we also show the best fits with the WDDl (green dashed 
line) and WDD2 (blue dotted line) SNIa yields, All abundances are 
relative to the Solar abundance standard of Asplund et al. (2004). 



fied by eye in the lightcurve from a low surface-brightness 
region of the CCDs and data from these intervals were 
excised. Point sources were detected in the 0.3-7.0 keV 
image with the wavdetect CIAO task. All sources were 
confirmed visually, and appropriate elliptical regions con- 
taining ^ 99% of the source photons were generated. 
Data from these regions were excluded in subsequent 
analysis. 

In Fig m (top right), we show a smoothed, flat- fielded 
Chandra image, having re moved the point so urces with 



we plot (data — model)^ / model , corresponding (approx- 
imately) to the residuals from this fit. To bring out 
the structure, we smoothed this image with a Gaussian 
kernel of width 3 pixels. The lack of significant, coherent 
residuals indicates that the X-ray image is very relaxed. 
An alternative view of the central (< 15") region is given 
in the lower left panel, which shows the azimuthal varia- 
tion in the surface brightness, integrated over two radial 
bins. There are no statistically significant residuals from 
an elliptical /3-model fit. Although the central isophotes 
in the smoothed Chandra image do appear slightly more 
flattened, this feature is not statistically significant, and 
may represent a statistical fluctuation. 

Spectra and associated flux-weighted responses were 
extracted in a series of contiguous, concentric annuli 
placed at the X-ray centroid. The widths of the annuli 
were chosen to contain approximately the same number 
of background-subtracted counts, while ensuring suffi- 
cient photons for useful spectral analysis. The result- 
ing annuli had widths larger than ~3", which is suffi- 
cient to prevent spectral mixing between adjacent annuli 
on account of the finite spatial resolution of the mir- 
rors. Data in the vicinity of point sources and chip gaps 
were excluded. The spectra were fitted simultaneously 
with Xspec, similarly to the XMM data, and the fit was 
restricted to the 0.5-7.0 keV band. To model the back- 
ground, we used a similar approach to our XMM analysis, 
although the surface brightness of the background com- 
ponents was assumed to be constant over the region of 
interest (which is entirely confined to the S3 chip). Rep- 
resentative spectra, and the best-fitting models, result- 
ing from a joint fit to all the Chandra or XMM data, are 
shown in Fig ^ The best-fitting abundance profile, and 
abundance ratios are given in Fig ^ while the projected 
temperature and density profiles are shown in Fig |[ 

4. MASS MODELLING 

Under the hydrostatic approximation, we transformed 
the projected density and temperature data into mass 
constraints with the e ntropy-based "for ward fitting " 

W ), 

errors. 



(2011 



technique described in Humphrey et al 
which enables tight control over systematic 
in comparison with other popular methods, such as 
"smoothed in yersion" (for a review of n iass modelling 



methods, see Buote fc Humphrey 2012a). Briefly, the 



entropy-based forward fitting method involves solving 
the equation of hydrostatic equilibrium to compute tem- 
perature and density profile models, given parametrized 
mass and entropy profiles. The models were then pro- 
jected onto the sky and fitted to the projected temper- 
ature and density profiles, taking into account covari- 
ance between the density data points. As is standard, 
we assumed spherical symmetry, which does not intro- 



^ , ^ duce substantial biases into the inferred mass distribu- 

the algorithm outlined in Fang et al. ( ^2009 ). The im- tion ( Buote fc Humphrey 2012c ). So as not to violate 



age was smoothed with a Gaussian kernel, the width of 
which varied with distance from the nominal X-ray cen- 
troid according to an arbitrary powerlaw, ranging from 
^1" at the centre of the image to ~0.9' at its edge. The 
image is smooth, albeit slightly elliptical. To search for 
more subtle structure we used dedicated software to fit 
an elliptical beta model (with constant ellipticity) to the 
unsmoothed (flat-fielded) image. In Fig|| (bottom right), 



the Schwarzschild criterion for stability against convec- 
tion, the entropy profile must rise monotonically, so we 
parametrized it as a constant plus a broken powerlaw 
model. We explored different models for the mass dis- 
tribution, as discussed below, and allowed the parame- 
ters describing the entropy profile, plus the logarithm of 
the gas density at a fiducial radius, and the mass pro- 
file parameters to fit freely. Parameter space exploration 
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Fig. 6. — Radial temperature (top panels) and density (bottom panels) profiles for NGC1521. We show the projected Chandra profiles 
in the left column (triangles), the projected XMM profiles in the right column (stars). The light blue region indicate the effect on the 
XMM density and temperature profiles of different data-analysis choices; the error-range shown combines both systematic and statistical 
uncertainties. Overlaid (solid lines) are the best hydrostatic model fits to each dataset, which match the data very well. The dotted lines 
are the best models if dark matter is omitted, for which the fit, in particular to the XMM data, is very poor. 

plying DM is required at 8.2-a. In Table |l|, we give the 
stellar M/L ratio, and show the marginalized total mass 
within i?A for various different overdensities (A), and 
the corresponding concentrations (ca — RA/fg, is the 
characteristic scale of the NFW model). 

In order to explo re how NGC 1521 sit s on th e a-Rc rela- 
tion discovered by Humphrey & Buotc (2010), we exper- 
imented with fitting the data within the central 50 kpc, 
under the assumption that the mass density profile is 
given by p oc r~". The model fits the data comparably 
to the NFW-|-stars model discussed above, confirming 
that the central part of the mass profile is approximately 
powerlaw in form. The constraints on the mass model are 
summarized in Table|^. We obtained a — 1.95tQQg, com- 
parable to other systems in Humphrey & Buotc ( 2010| ) 
with K-band Ro similar to NUC 1521 (=4.7 kpc, as given 
in 2MASS)0. 

5. SYSTEMATIC ERROR BUDGET 

In this section, we address the sensitivity of our re- 
sults to various data analysis choices that were made. 
Since it is generally impractical to express these assump- 
tions through an additional model parameter over which 
one can marginalize, we adopted the pragmatic approach 
of exploring how our results changed if the assumptions 

1" We used the 2MASS measurement of Rc for NGC 1521 

n nnahlo a r'r.naiatprjt rnmparlann with the measurements frOm 

Humphrey fc Buotc (2010). 



emploved version 2.7 of the MultiNest Bayesian code^ 
( Feroz et al. 2009] ). Initially we adopted flat priors for 
each fit parameter, but subsequently explored the impact 
of this choice on the fit results (§ |^). 

For the mass models, we first considered the case of 
no DM. We included a stellar light component, based on 
spherically averaging the deprojected, triaxial r nodel for 
the I -band l ight di scussed in § 2.2. As shown in Buotc fc 
Humphrey (2G12c), the mass inferred from spherical hy- 
drostatic methods should be very close to the spherically 
averaged true mass. The mass-to-light (M/L) ratio was 
allowed to fit freely. We also included a (fixed) black hole 
with mass Mr p=3 x lO^Mg. c onsist ent with the Mbh- 



(T* relation of [Gultckin et al.| (P009| ) , give n the central 
stella r velocity dispersion a,: = 224km s^i ( |Faber et aL| 
1989). The fit to the density and temperature profiles 



was poor (Y^/dof=1 07/21), as shown in Fig ^. Next, 
we added an NFW ( [Navarro et"aLl|l997D DM halo. We 
allowed ^ogioMvir and /ogioCvir (the DM halo concentra- 
tion) to fit freely (restricting 11 < logi^Mmr ^ 15 and 
< logioCvir < 2), and obtained a formally acceptable 
fit (x^/dof=29.5/19; Fig ^). The measured radial mass 
distribution is shown in Fig |7[ 

The improvement to the fit when dark matter was in- 
cluded was highly significant; the ratio of the Bayesian 
"evidence" returned for the two cases is 1.8 x 10"^^, im- 

® http:/ /www. mrao.cam.ac.uk/software/multinest/ 
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TABLE 1 



Mass rrsiit.ts and rrror RTinriFiT 


Test 


M,/Li 
MqLq-i 


log M2500 
[Mq] 


log C2500 


log M500 
[Me] 


log C500 


log Mvir 


log Cvir 


Marginalized 
Best-fit 


2 4Q+0 26 
^■^y_0.42 

(2.55) 


12.40 ±0.08 
(12.40) 


(0.58) 


12.59 ±0.11 
(12.59) 


(0.88) 


1 9 70+0.15 

(12.76) 


1 1Q+0-14 

j^-j^y_o.i7 
(1.16) 



AAC 
AStars 
AH.E. 



+0 

ABackground ±0.13 
ASWCX 
AInstrument 



0.652 (±0.37) +0.008 (±0.09) -0.086 
-0.045 (±0.08) 
+0.16 (±0.09) 



/+0.21 

\^-0.28 

0.18 f+°Ai 



on I +0.37^ 
^•^ I -0.53 ) 



) +0.003 -0.073 ( 



+0.35 
-0.50 



(±0.10) 



0.095 (±0.42) 
0.22 (±0.31) 



-0.05 

-0.090 
h0.08 



AFit radius +0.06 (±0.27) +0.20 



+0.12 

-0.07 
-0.16 
-0.11 
-0.10 
-0.13 



A3d 

AFit priors 
ASpectral 
AEntropy 
AWeighting 
ADistance 
ACovariance +0 



.353 (±0.53) 
.39 (±0.21) 

f+0.24 
1^-0.33 
f+0.23 
0.45 



-0 
+0 
-0 
+0 

-0.158 (±0.44) ±0 (±0.09) 



194 
.03 



-0.079 (trii) 

+0.06 (±0.08) 
-0.043 (±0.06) 



V- 



0.006 



59 f+0.36\ 
45 1^-0.49^ 

1^-0.82 1 



.16 



-0.054 {trol] 

+0.05 /'+0.13^ 
-0.02 \^-0.09 J 



^^-0.23 
+0.02 (±0.20) 

+0.08 f+0.2l\ 
-0.12 V-0.29^ 

+0.11 ip,]^ 

-0.243 (±0.28) 
-0.232 (±0.20) 

+0.03 {tin) 
-0.296 (±0.14) 
+0.04 (±0.15) 
-0.010 (±0.19) 
+0.04 

±0.08 (tr^i 

+0.07 f+0.23^ 
-0.18 1-0.32 



1+0.003 

-0.043 (±0.11) 
+0.15 (±0.12) 

tro°7 (±0.14) 
-0.122 

+0.15 (±0.22) 
+0.24 (±0.17) 
-0.096 



0.073 [t.°ofr 
+0.18 

+0.00 f+0.21 

-0.03 \^-0.15 

(±0.21) 

+0.10 



-0.02 



/+0.20 



[-0.13) 

+0.15 (±0.13) 
-0.051 (±0.09) 
-0 018 I +0-15 



+0.16' 
0.10 
0.16'' 
0.08 , 
+0.09 f+0.20^ 
0.05 V-0.12^ 



-0.025 
-0.080 



-0.219 (±0.25) 

-0.212 (±0.18) 

+0.01 (1°;??) 
-0.275 (±0.13) 
+0.03 (±0.14) 

-0.010 (±0.17) 

+0.005 (±0.19) 

±0.08 

+0.07 f +0.22\ 
-0.16 1-0.31 J 



\-0.12 
±0.03 (±0.13) 

+0.17 (±0.13) 
tao3 (±0.16) 

-0.104 

■ f0.30'\ 
-0.21 j 



-0.072 



+0.18 
-0.26 



+0.18 



+0.00 

-0.04 



+0.21 
-0.14 



tnt (±0.20) 



^0.18 



hO.lO 
-0.206 



+0.16 
-0.21 
/+0.21 
1-0.26 



+0.28 (±0.19) 
-0.084 



f+0.26^ 
I -0.09 , 



-0.19 (±0.14) 
-0.027 (±0.10) 

+0.01 



0.198 (±0.17) 
/+0.3A 

(-0.191 



+0.03 



-0.015 

+0.00 /+0.19" 

-0.06 l-O.ll 

+0.13 /+0.24 
-0.04 1-0.12 



-0.261 (±0.13) 

+0.03 (±0.13) 
-0.028 (±0.17) 
+0.03 (±0.19) 
±0.08 {till) 

+0.06 f+0.2l\ 
-0.15 1-0.30 ) 



Note. — Marginalized values and l-cr confidence regions for the stellar mass-to-light (M*/Lk) ratio and the enclosed mass and 
concentration measured at various overdensities. Since the best-fitting parameters need not be identical to the marginalized values, 
we also list the best-fitting values for each parameter (in parentheses). In addition to the statistical errors, we also show estimates 
of the error budget frorti- possible sources of systematic uncertainty. We consider a range of different systematic effects, which are 
described in detail in § pi specifically we evaluate the effect of the choice of dark matter halo model (ADM), adiabatic contraction 
(AAC), treatment of the stellar light (AStars), plausible deviations from hydrostatic equilibrium (AH.E.), treatment of the background 
(ABackground) and the Solar wind charge exchange X-ray component (ASWCX), the instrumental inter-calibration (AInstrument), 
the radial coverage of the data being fitted (AFit radius), deprojection (A3d), priors on the model parameters (AFit priors), spectral 
fitting choices (ASpectral), the parameterization of the entropy model (AEntropy), removing the emissivity correction (AWeighting), 
distance uncertainties (ADistance), and the treatment of covariance between the temperature and density data-points (ACovariance). 
We list the change in the marginalized value of each parameter for every test and, in parentheses, the statistical uncertainty on the 
parameter determined from the test. Note that the systematic error estimates should not in general be added in quadrature with the 
statistical error. Since there is no theoretical interest in the distribution of scale radii for the cored logarithmic DM model, we do not 
include changes on the concentration in our error budget for this choice of profile (ADM profile). Likewise, the tests involving different 
stellar light modelling approaches (AStars) involved using light profiles from different optical filter, and so we omit the M.^/L/ ratio 
from the error budget for that choice. 



were adjusted in an arbitrary, but representative, way. 
We focused on those systematic effects likely to have the 
greatest impact on our conclusions. In Tables we list 
the change in the marginalized value of each parameter. 
We discuss each test in more detail below. In summary, 
most of the inferred systematic errors are comparable to 
the s tatistical errors. For a couple of the tests, marginally 



the adopted priors (a flat prior on the asymptotic circular 
velocity, between 10 and 2000 km s~^, and a flat prior 
on logioVc^ where Tc is the core radius, over the range 
< logioVc < 3.) is a poorer description of the data at 
-3.4-cr. 

A theoretical modification to the DM p rofile is ex- 



significant (~2-cr 
reductions m tb.v 



increases in Mvn, 
are seen (by as much as 
our conclusions are not strongly affected. 



and corresponding 
U.UVj, but 



pected to arise from adiabatic contraction j Blum enthal 



et alj|l986t [Gnedin et alj|2007| |Abadi et al.||201(j| ), which 
causes the halo to become cuspier due to the gravita- 



5.1. Dark matter halo 

The accurate computation of the gas distribution out 
to Rvir is contingent upon the accurate modelling of the 
gravitating mass distribution. While the NFW DM halo 
model is well-motivated theoretically, we also considered 
the "cored logarithmic" model that is sometimes used iii 



stellar dynamical studie s (e.g. Binney & Tremaine 
fien fc Gebharf^ pO10| ). This model requires a 
higher Mvir, and correspondingly lower fb,vir ( "ADM pro 



2008 
ightly 



file" in Tables |] and H). Although we cannot distin- 
guish between the NFW and cored logarithmic model on 
the basis of alone, the ratio of the Bayesian evidence 
(7 X lO^**) implies that the cored logarithmic model, with 



tional influence of the bary ons. Modifyi n g the NFW 
proflle with the algorithm of pnedin et al. ( 2004| )P^, has 
only a very slight effect on the best-fitting mass model 
("AAC" in Tables |l| and |), except for the stellar M/L 
ratio, which is reduced to make room for the increased 
DM fraction predicted in the inner parts of the galaxy. 

5.2. Stellar light 

Accurately decomposing the gravitating mass distri- 
bution into the luminous and dark components requires 
an accurate model for the (deprojected) stellar light. In 
the case of NGC1521, we used a triaxial deprojection 

Using the CONTRA code publicly available from 
http: / /www. astro. lsa.umich.edu/ ~ognedin/contra/ 
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B-band surface brightne ss profiles along the ma jor and 



minor axis published by Capaccioli et al. (1988) with a 



Fig. 7.— Radial mass profile of NGC 1521. The solid (black) 
line indicates the total enclosed mass (and the grey shaded region 
indicates the l-cr error in the total mass distribution), the dashed 
(red) line indicates the stellar mass, the dotted (blue) line is the 
dark matter, and the dash-dot (magenta) line is the gas mass con- 
tribution. Overlaid are a set of data-points drriirrd from n mnrf 

<:raHitiririal "smnnthpri iyrerginn" approach ([Buote fc Humphrey 

p012al; Humphrey et al. 2009a). We stress that the model is not 
nttca to tnesc aata, out is aerived independently from the temper- 
ature and density data. 



TABLE 2 

Mass slope results and error budget 



Test 


log M75 


a 




[Mq] 





AH.E. 

A Background 
ASWCX 
AInstrument 
AFit radius 
A3d 

AFit priors 

ASpectral 

AEntropy 

AWeighting 

ADistance 

ACovariancc 



-1-0.03 (±0.05) 
laoa (±0.06) 
+0.01 {tl°l) 
t°o:ll (±0.06) 
+0.03 (iHf 
-1-0.05 (±0.07) 
-o°ol (±0.06) 
-0.020 (±0.04) 
-0.023 (±0.05) 

0.03 flH^) 



+0.04 (±0.05) 
to;o2 (±0-05) 
-0.010 (±0.05) 
t°ofo (±0.06) 



-0.028 



f +0.03 



-0.037 

+0.05 



h0.05 
-0.06 



(±0.07) 



\^-0.05 

-0.075 (±0.06 

t°o°ol (±0.06) 
+0.008 (±0.04) 
+0.02 (±0.05) 

-0.033 (IHI 
+0.03 (±0.05) 

±0.04 



Note. — Marginalized values and l-cr confi- 
dence regions for the enclosed mass at 75 kpc 
(.^^^75) and the negative logarithmic slope of the 
mass profile (a), when fitting only a single pow- 
erlaw to the total mass distribution. We also pro- 
vide the best-fitting parameters in parentheses, 
and a breakdown of possible sources of systematic 
uncertainty, following Table hi 



procedure (§ 2.2). Although we averaged this profile 
spherically for use with our X-ray modelling code, we 
have previously shown that this ap proach does not gen- 
erally introduce significant biases ( Buote fc Humphrey 
^012b| Jc|). Nevertheless, the deprojection procedure may 
not be unique, and so it is important to explore how 
sensitive our results are to the details of the stellar light 
modelling. 

We experimented with three alternative prescriptions 
for deprojecting the stellar light. First, we fitted the 



model comprising 6 multiple, concentric 3-dimensional 
Gaussian density distributions, that were projected onto 
the sky, assuming an edge-on, oblate geometry. This ig- 
nores the isophotal twist. Second, we fitted the K-band 
2MASS image in the central ~4' region with a model 
comprising an elliptical S ersic model, which wedepro- 
jected with the formula of Prugnicl & Simie ^(|199^), as- 
suming an edge-on oblate spheroidal geometry, i'lnally, 
we adopted a spherical, deprojected de Vaucouleurs 
model, the effective radius and luminosity of which were 
set to match the catalogued 2MASS K-band values. We 
found that the treatment of the stellar light primarily 
only affected the total stellar mass, so that fb,2500 is very 
sensitive to this choice. However, by Rvh- the impact is 
much less significant ("AStars" in Tables ^ and H). 

5.3. Hydrostatic equilibrium 

Although the gas in morphologically relaxed early-type 
galaxies is expected to be close to hydrostatic, recent 
work suggests that a small amount (<30%) of nonther- 
mal support is present in the very central parts of some 



et 



ail ^008| ; |Das et al.| |201(]| ; 



galaxies fe.g. Sliurazoy 

Humphrey et aTT 2012b lb investigate whether devi- 
ations from hydrostatic equilibrium at this level would 
quantitatively affect our conclusions, we modified the hy- 
drostatic equation used in our modelling code to include 
a plausible nonthermal component. First we considered 
a nonthermal pressure fraction profile similar to that in- 
ferred for the galaxy NGC 4649, which was fixed at ^25% 
at the centre and fell to ^10% by 20 kpc, vanishing out- 
side 30 kpc (Humphrey ct al. p012b ). This is similar to 
the implied nonthermal pressure profiles that have been 



inferred in (albeit a handful of) other systems (e.g. Das 



t al. 2010). Since the nonthermal pressure is most im 
portant in the central few kpc, we found that adding 
this component did not significantly affect the global pa- 
rameters of the system, although the stellar M/L ratio 
was increased by ~18% to compensate. As a more ex- 
treme alternative, we also considered a uniform nonther- 
mal pressure fraction fixed at 25%, similar in magnitude 
to that inferred fr om previous studies o f the central parts 
of of galaxies (e.g. [ Churazov et al. |200S , [2010 ; Humphrey 
^t al.|^0094 |2012b| ; pas et alj 12010^ In this case, the 
stellar M/L ratio was mcreasea more substantially (by 
~33%), and the global mass raised by ~0.15 dex. Con- 
versely, the baryon and gas fractions at fixed overdensity 
were found to be relatively insensitive to this choice, since 
both stellar and gas mass are increased, while the slightly 
larger Rsqd, for example, encloses more gas ("AH.E." in 
Tables |l|-|). 

5.4. Background 

Since the data were background-dominated in the 
outer XMM annuli, the treatment of the background was 
a potentially serious source of systematic uncertainty. To 
investigate the extent to which our results are sensitive 
to this, we explored a number of different choices in our 
background treatment. First, for the Chandra data, we 
adopted the standard blank-field events files distributed 
with the CALDB to extract a background spectrum for 
each annulus. Since the blank-field files for each CCD 
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TABLE 3 

Baryon fraction results and error budget 



Test 


/9,2500 


fg, 500 


fg,vir 


/6,2500 


/t.soo 


fb.vir 


Marginalized 
Best-fit 


0.019 ± 0.002 
(0.018) 


(0.039) 


0.10 ±0.02 
(0.095) 


0.16 ± 0.03 
(0.178) 


0.13 ± 0.03 
(0.142) 


0.16 ± 0.03 
(0.165) 



ADM profile 
AAC 
AStars 
AH.E. 

ABackground 
ASWCX 
AInstrument 
AFit radius 
A3d 

AFit priors 

ASpectral 

AEntropy 

AWeigliting 

ADistance 

ACovariance 



±0 (±0 
±0 (±0 
±0 (±0 
-0.003 
-0.007 
+0.001 

±0 (1° 
-0.004 
+0.004 
±0 (±0 
+0.002 



-0.007 (±0.005) 



.002) 
.002) 

.003) ±0 (±0.01) 
(±0.002) 
(±0.003) 
(±0.003) +0.004 (±0.01) 

.oo2^ 

.004 J 
/'+0.003^ 
1^-0.003 J 



-0.001 (tlZ] 



-0.007 (±0.01) 
tlTs (±0.01) 



-0.006 



(±0.006) 

.003) 

(±0.002) 



±0 



/+o. 



V-o 
-0.003 

±0.003 

±0 (±0 



/'+0.003^ 
1^-0.002 ) 

/'+O.OO3S 
1^-0.002 j 

,002) 



-0.011 (t^Z 

+0.02 {tVol) 
-0.005 (±0.01) 
+0.004 {t°o fol) 
±0 (±0.01) 
-0.006 (±0.01) 



+0.004 
-0.006 
+0.002 
-0.004 



+0.01 

-0.01 
+0.01 
-0.01 



-0.029 (±0.01) 
-0.007 (±0.02) 
+0.009 (±0.02) 
-0.015 (±0.02) 
tVol (±0.03) 
+0.007 (±0.03) 
-0.023 (tH^ 
-0.059 (t^ot 
+0.05 

-0.019 (±0.02) 
+0.01 {tlf, 

'+0.02 
-0.03 
+0.02 
-0.02 



+0.005 (±0.04) 
-0.042 (±0.02) 
-0.085 



-0.003 
-0.012 
±0.01 (±0.02) 

+0.005 /'+0.03^ 
-0.012 I -0.03 J 



±0.02 

+0.01 

-0.04 

+0.006 
-0.029 
-0.056 

-0.036 

+0.00 
-0.01 

±0 
-0.013 
-0.019 
±0.03 
-0.017 



+0.03 
-0.02 

(±0.03) 
(±0.04) 

/+0.04 
1^-0.04 

(±0.03) 

+0.03 
-0.02 

+0.05 
-0.03 

(±0.03) 
(±0.02) 

-0.03 
-0.02 

(±0.03) 
(±0.03) 
(±0.03) 



-0.021 (±0.02) 
-0.026 (tH^ 
-0.049 (±0.02) 
(±0-03) 
1^4 (±0-03) 
+0.009 flHa) 



-0.025 (iHi) 
-0.049 (±0.03) 
+0.002 (±0.05) 
(±0.03) 
+0.01 (tlf,) 
-0.005 (±0.03) 
-0.012 (±0.03) - 

to.02 (±0.03) 
-0.010 (±0.03) -0.015 (±0.04) 



-0.056 [tVol 
-0.023 (±0.03) 
-0.029 (±0.03) 
iao2 (±0.04) 
laoe (±0-05) 
+0.02 (±0.05) 

-0.046 (±0.05) 
-0.071 (±0.04) 

+0.04 (ini) 

-0.030 (±0.04) 
+0.02 (tl f, 
-0.009 (±0.03) 
-0.015 (iHi 

+0.03 /+0.04 
-0.01 1-0.04 



Note. — Marginalized values and 1-cr confidence regions for the gas fraction (/g,A) and baryon fraction (ft^A) measured 
at various ovcrdensities (A). We also provide tha_best-fitting parameters in parentheses, and a breakdown of possible 
sources of systematic uncertainty, following Table hi 



have different exposures, spectra were accumulated for 
each CCD individually, scaled to a common exposure 
time and then added. The spectra were renormalized to 
match the observed count-rate in the 9-12 keV band. 
These "template" spectra were then used as a back- 
ground in Xspec, and the background model components 
were omitted from our fit. This did not strongly affect 
our conclusions. 

Since the non X-ray component dominates the XMM 
background for much of the band-pass, we explored an 
alternative means of accounti ng for it. Specifically, we 



adopted the ESAS a lgorithm (Kuntz & Snowden 



2008 



^nowdcn et al. 2008 ). This involves choosing non X-ray 
background template files that match the count-rate and 
hardness ratios of photons in the unexposed portions of 
the CCDs. We added the X-ray background and instru- 
mental lines, as in our standard modelling procedure. 
Using this approach, we found a modest reduction in the 
gas density at all radii (Fig In the innermost regions 
this arises due to a slightly higher best-fitting abundance 
(Zfc~ 0.6). Although this lowers the baryon fraction 
slightly (by 4% at R500), this change is comparable to the 
statistical error. Next, we explored the impact of varying 
the X-ray background. We varied the slope of the cosmic 
X-ray component (powerlaw model) by ±5%, changed 
kT of the local hot bubble component to 0.1 keV, and 
varied kT for the galactic hot gas component between 
0.18 keV and 0.24 keV, which span the range of tem- 
peratures obtained from fitting nearby blank-sky fields 
(§ 3.1). Although these choices had a measurable impact 
on the global parameters we derived for NGC 1521, our 
conclusions were largely unaffected. We summarize the 



results in Tables 0-| ( "ABackground" ) . 

5.4.1. Solar Wind Charge Exchange 

A potentially relevant, time variable, soft background 
component can arise due to the interaction of the So- 
lar wind with interstellar material and the Earth's exo- 
sphere. To explore the possible importance of this "Solar 
Wind Charge Exchange" (SWCX) compo nent, we used 



data from the WIND-SWE experiment (Ogilvie et al 
1995)^ to ide ntify periods of s trong Solar wind activ- 



(2004), we assumed the 



Tty^ Following Snowden et al. 
SWCX component is negligible for a Solar wind proton 
flux level measured to be < 3 x lO^cm^-^ s~^, and excised 
all other data. While this eliminated ^-^30% of both the 
Chandra and XMM exposures, we found that this choice 
only had a minimal inipact on the derived parameters 

indicating that the SWCX 

em in our analysis. 



("ASWCX" in Tables _ 
is not a significant prob. 



5.5. Instrumental inter- calibration 

In our default XMM analysis, we considered only the 
MOSl and M0S2 instruments, as we have found the 
calibration of the PN to be more uncertain in the low- 
temperature (kT'^0.5 keV) regime. Nevertheless, given 
its large collecting area, it is important to investigate 
whether the PN data can add information in the low sur- 
face brightness outer regions of the galaxy. To explore 
this, we first fitted the PN spectra, similarly to the MOS 
data, to obtain the gas temperature and density profile 
out to ^--^200 kpc. We then fitted these profiles in tan- 
dem with the Chandra data. Although the derived Mvir 



http; / /web.mit.edu/spacc/virwvir/virind_data.html#Protons 
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was slightly larger (and, consequently, fgas and fb val- 
ues were smaller) than our best fitting case, the results 
were mostly consistent within errors ( " AInstrument" in 
Tables 

5.6. Radial coverage 

The constraints on the global galaxy properties, es- 
pecially those derived from a full mass decomposition, 
can be sensitive to th e radial range being fitted (e.g. 
[Gastaldello et al. 2007|) . To explore the sensitivity of our 
results to this choice, we restricted the radial range in 
our fit to < 100 kpc (by excluding the outer three XMM 
datapoints). This increased the mass of the system by 
0.28 dex (simultaneously reducing fgas and fb), but only 
at the expense of larger error bars ( " AFit radius" in Ta- 
bles and I). 

For the powerlaw fits to the total mass distribution, by 
default we excluded all data outside '^50 kpc. Since the 
mass profile must deviate from a pure powerlaw at large 
scales, we explored the impact on our results of extending 
to slightly larger scales in the fit (out to ~70 kpc). This 
did not affect our results significantly ("AFit radius" in 
Table I). 

5.7. Projection/ Deprojection 

In this work, we fitted the projected, rather than the 
deprojected data. To determine the mass profile, we 
modelled the projected temperature and density in each 
annulus by evaluating the hydrostatic model for the tem- 
perature and density in three dimensions, and projecting 
it onto the line of sight. In general, this procedure leads 
to smaller statistical error bars, but may introduce ad- 
ditional systematic uncertainties relat ed to how the pro- 
iected quantities are computed (e.g. pastaldello et aL 
^007 ). To explore the likely impact of these effects, wc 



performed a spherical deprojection by adding multiple 
"vapec" plasma models in each annulus, w ith the rel 



ative normalizations tied appropriately (e.g. Kriss et al 
1983 ). While functionally identical to the "projct" model 



in Xspec, this allowed data from both MOS detectors to 
be fitted simultaneously. Given the amplific ation of noise 
by the deprojection procedure (e.g. § 3.3 of Buote 2000a ; 
Finoguenov fc Ponman| 1999 ), we excluded the (noisy) 



outer three XMM annuli from this analysis. To account 
for emission projected into the line of sight from regions 
outside the outermost annuli, we added an APEC plasma 
component to the spectral model with an abundance 0.2 
(consistent with the outermost annuli) and the tempera- 
ture and normalization determined in each annulus from 
projecting onto the line of sight the best-fitting gas tem- 
perature and density models evaluated beyond the outer 
bin. 

Although we found no evidence of biases when com- 
paring the projected and deprojected results, the reduc- 
tion in radial range fitted, in conjunction with the am- 
plification of error bars during deprojection, resulted in 
much poorer constraints on the global quantities com- 
puted within Rvir- Nevertheless, at smaller scales (e.g. 
R.2500): the results were much more robust ("A3d" in 
Tables |-|). 

5.8. Priors 

Since the choice of priors on the various parameters is 
arbitrary in our analysis, it is important to determine to 



what extent they could affect our conclusions. To do this, 
we replaced each arbitrary choice in turn with an alter- 
native, reasonable prior. Specifically, for each parameter 
describing the entropy profile, we switched from a fiat 
prior on that parameter to a flat prior on its logarithm. 
We used a flat prior on the DM halo mass, rather than on 
its logarithm, and, instead of the flat prior on log cdm, 
we adopted t he di stribution of c around M found by 
Buote et al. (2007) as a (Gaussian) prior. We also re- 
placed the flat prior on the M*/i/ ratio with a Gaussian 
prior, corresponding to the best-fltting M/L ratio and 
err or-b ar derived from fitting the published Lick indices 
(§ |6j| ) . For the single powerlaw mass model fits (Ta- 
ble^), we used a flat prior on M75, rather than its loga- 
rithm, and a flat prior on the logarithm of a. The effect 
of these choices is comparable to the statistical errors on 
each derived parameter ("AFit priors" in Tables 1 and 

i- 

5.9. Other tests 

We here outline the remaining tests we carried out, as 
summarized in Tables |l| and First of all, to assess the 
sensitivity of our results to the choice of plasma code em- 
ployed, we experimented by replacing the APEC plasma 
model with a MEKAL model. This had very little impact 
on our conclusions ("ASpcctral" in Tables 0-||). 

In order to explore the sensitivity of our results to 
the entropy parameterization we adopted, we experi- 
mented with allowing an additional break at large ra- 
dius. We found that the break radius was poorly con- 
strained (35^23 kpc) and adding it did not significantly 
improve the fit (Ax^=3 for 2 d.o.f.). Based on the ratio 
of the Bayesian evidence (6 x 10^^), the model without 
the break is actually preferred at the level. In- 

cluding the break had a minimal effect on the best-fitting 
derived parameters ("AEntropy" in Tables 

In our default analysis, the projected temperature and 
density profile were weighted by the gas emissivity, folded 
through the instrumental response s (for details, see Ap- 
pendix B of Gastaldello et al. 2007). Since the computa- 
tion of the gas emissivity assumes that the three dimen- 
sional gas abundance profile is identical to the projected 
profile (which is unlikely to be true), we explored the sen- 
sitivity of our results to this approximation by adopting 
the extreme approach of ignoring the spatial variation of 
the gas emissivity altogether. We found that this had a 
very small effect on our results (AWeighting in Tables 111- 

To examine the error associated with distance uncer- 
tainties, we varied the distance to NGC 1521 by ±20%, 
finding the effect, particularly on fb and fgas to be rela- 
tively minor ("ADistance"). Finally, to examine the pos- 
sible errors associated with our treatment of the covari- 
ance between the density data-points, we investigated 
adopting a more complete treatment that considers the 
covariance between all the temperature and density data- 
points, as well as adopting the more standard (but incor- 
rect) approach of ignoring the covariance altogether. We 
found that this had a non-negligible, but modest impact 
on the derived parameters ( " ACovariance" in Tables n[~ 
i. ^ 

6. DISCUSSION 
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6.1. The EllXr Galaxy Survey 



pait 



to pro- 
hot gas 
axy-scale 



The EllXr Galaxy Survey is intended, in 
vide a sample of candidate L* galaxies witl 
halos which are likely to be hosted in a ga ; 
f<10"Mo) DM halo. The full sample ( piT 
^012 ) includes several objects that are known 
these characteristic s fe.g. NGC 720: jHum phj-cv et al 
[Humphrey et al" [jUU6 



)tc et al 



2011; NGC4125: 



[J 'Sullivan et al 



et al 



pOOT] ), but the majority of 
ies have not been studied before in the X-ray. 



GC 7706: 

galax- 
1KGC1521 



was identified aa an X ray bright object, baa sd on the 
shori initial EllXr XMM observation , and tar jetted for 



deep follow-up with both Chandra and XMM. 

The X-ray morphology is very relaxed, withoi t 
disturbance and there is no bright AGN fas evir 



modest radio emission: Condon et al. 1998, an 



of a central X-ray point source), suggesting nc 
of jet / gas heating in the recent past. The hot 
clearly detectable out to ~200 kpc (-80% of R, 
is 240 ± 22 kpc), and the temperature profile 
imately isothermal at kT~0.5 keV. Despite tl 
higher mass, the properties of NGC 1521 are re 
similar to NGC 720, making NGC 1521 one of 1 
est DM halos for which interesting X-ray cons 
the mass profile have been found, and helping 
late the very sparse low-mass end of the c 
with vital new data. The properties of NGC 
in fact, very close to our expectations for ElIXi 
providing a crucial validation of our observing 



6.2. 



Hydrostatic Equilibrium 

The best-fitting hydrostatic model fits the d« nsity and 
temperature data-points well, consistent witli the gas 
being close to hydrostatic. Despite nontrivial temper- 
ature and density profiles that cannot be paraneterized 
by simple models individually, a smooth, physical mass 
model, coupled to a monotonically rising entropy profile 
(i.e. that is stable against convection), was able to repro- 
duce them well. If the gas is far from hydros ;atic, this 
would require a remarkable conspiracy betweer the tem- 
perature, density and inferred mass profiles, ^^he close- 
ness of the system to hydrostatic is unsurprif ing given 
its relaxed X-ray morphology (§ ||). 

Nevertheless, modest deviations from hydrosi atic equi- 
librium cannot be entirely ruled out. While 
structure formation simulations suggest that 
around galaxies sh ould, indeed, be quasi-hydrof 



prain et al. 2010), deviations from the hydro 



proximation may introduce systematic errors on the re- 
covered mass of as much as —25%, if they a 'e similar 



gai et alj p007| ; [Piffarctti fc Valdarnini p008 
2009 ). In order to constrain such non-nydrot 



et al 



Eit 



obvious 
ced by its 
i the lack 
evidence 
ias halo is 
,-,00, which 
approx- 
e slightly 
markably 
he small- 
iraints on 
to popu- 
vir plane 
1521 are, 
galaxies, 
strategy. 



Jl 



lumerical 
hot halos 
tatic (e.g. 
tatic ap- 



to clusters (e j :. Jrsai et al.| |l994 iRasia et al.| 2006|: Wa- 



ang et al 



observationally, recent studies have begun to 



ydrosta tic ett'ects 



the mass distributions of galaxies inferred from stellar dy 
namics modelling to those independently obtai 



rays, including several systems t hat are manifestly more 



disturbed than NGC 1521 f e.g. Humphrey ct 
b008l 



Churazov et al. 



3hen fc Gebhardtll2010 



2010; Gcbhardt k Thomas 



al.| t20l2b| ). Altfen^ 



al 



0as et al.||201^, |011 
dI o 



2008 



2009 



1 better control of systc 
rors may be necessary before deviations from h; 
equilibrium can be measured accurately (e.g 



umphrey 



matic er- 
"drostatic 



Churazov 



red by X- ( |200'i ) reported substantially enhanced concentrations 



Hum jhrey et ai 



unde: 
tral r 
mal 
and 
mass 
is prdsent 



])ressure profiles consistent with these observations, 
fpund that, aside from a modest increase in the total 
of the system, if nonthermal pressure at this level 
in NGC 1521, our conclusions about the global 
prop(hrties of the system were largely unaffected 

While subjec t to uncerta i nties in the stellar initial mass 
func t 
roy 



2010) and the star-formation history of the galaxy. 



complarisons between the stellar M/L ratio inferred from 
e stellar population (SSP) models and from X-ray 
E!S similarly suggest that de viations from hydrostati c 
briu m are generally sniall fjHumphrey ct al. 2009a). 



aje 



simp 
studi 
equil 
Folio 
lar 

NGO 

to 

lishe( I 
and 
to ar 



vmg 



the 



agree 
This 
supp 
good 
ties 



by lOgando et al] f|2008|). For a|Kroupa| ( |200l| ) IMF 
SSP model of [Maraston (2005), this corresponds 
I-band stellar M/L ratio of 2.6±0.6MqLq\ which 

) well with our measurement (2.491o'42M0Lq"'^). 
s consistent with observations of other galaxies, and 
jrts the idea that the hydrostatic approximation is 
for NGC 1521. Nevertheless, given the uncertain- 
( Doth statistical and systematic) in this comparison, 
nontltcrmal pressure providing as much as —40% of the 
support could be consistent with the dat a, partic 



total 
ularb 
for a 
is wi 
precife' 
will 1 



XMM 
matt' 



and 
ray 



(1200^ 1) 



at IC 
mass 
more 
is e 
tions 



Mace 



m a 
Whil 
theri 
nor 
ture 



(200' 



2008 



Gebhardt fc Thomaiij |2009| ; pas et al.||20Tl| ; 
2012b), the X-ray measurements may 
estimate the mass by as much as —30% in the cen- 



5-10 kpc. In § 5.3, we explored plausible nonther- 



(e.g. Treu et al. 2010; [van Dokkum fc Con- 



( poof 



Humphrey et al 
(8 . S^i g Gyrj and metallicit y ([Z/H]=0.30; !^ 



we estimated the stcl- 



1521 by fitting the models of [Thomas et al 
H /3, Fe5270, Fe5335 a nd 



Ml 



^ ( |2003| ) 
b Lick indices pub 



if the IMF is more bottom heavy than Kroupa; 
Salpeter IMF, for example, M/L=4.0 ± 0.9, which 
ihin — 2-tT of our measurement. To quantify more 
>ely the accuracy of the hydrostatic approximation 
kely require sophisticated (orbit-based) stellar dy- 
;al modelling of this triaxial galaxy. 

6.3. Mass profile 

Ba sed on our hydrostatic analysis of the Chandra and 
data, we were able to confirm the presence of dark 
r in NGC 1521 at high significance (8.2-cr). The 
constlraints our model was able to place on the virial mass 
( onccntration (Fig |^) are competitive with other X- 
c etermined measu rements in this mass range (e.g. 
Humphrey et al. 2006). The system i s clea rly consistent 
both with the empirical Buote et al. ( ^007) Cy jr-M 
latioi 



and also the theoretical model of 



Maccio et al 



(which predicts a slightly lower logiocvir — 0.97 
^■^Mq). This may slightly ease tension at the low 
end between theory and observations, although 
observations are needed to resolve this issue. It 
:xj|)ected that isolated systems have higher concentra- 
on average, than the population as a whole, al- 
compare though the effect is not predicted t o be dramatic fe.g 



io et al. 2008 ). Nevertheless, Khosroshahi et al 



small sample of fossil groups and isolated galaxies. 
NGC 720 could be consistent with this picture, 
is no evidence of such an effect in NGC 1521, and 
d oes the EllXr ga l axy N GC 4125 show such a fea- 



Humphrcy ct al. 2006). A simil ar conclusion was 



react ed for the fossil groups studied by pastaldello et al 
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Fig. 8. — Concentration-virial mass constraints for NGC1521, 



shown ainn 



Buotc ct al 
nonrnjirs lor 



with. 



(2011) and Humphrey et al 



rn 


iphrpv 


( 





the sample of galaxies, groups and clusters in 

(2007). Where pnaaihlp wp ahnw the 1-^ rnnfiHenrr 

the nhjerts in [Hi|n:)phrP3[ et al (2006), Humphrey et al. 

~ I hp snlirl hhip line is thp 

Buote et al (2007), and the 
blue shaded region indicates the l-cr intrmsic scatter ot the points 
about it. NGC 1521 is in good agreement with this relation. 




Fig. 9. — Best-fitting (total) mass slope (a) versus K-band Rg re 

lafinri for NCiC 1 .^^91 (hlnp gt.ar'l Wp mfprlay (hlapV pirelpa'l thp Hat.a 

from,fJiLmphrey & Buotc (|2010), and the data derived by Auger 
et al. 




|J trom strong lensmg 



and stellar kinematics within ~K,c 



(grey stars) . The snliH and HnttpH (rprl^ linpa are the predictions 
of the toy model of Humphrey & Buotc (2010) for different Sersic 



indices (n) for the stellar ugnt 

NGC 1521 is a more luminous galaxy than NGC 720 
(by ~60% in the K-band, and ~90% in B), and so it 
is not surprising that it sits in a more massive halo (by 
a factor ^^2). The allowed mass range of NGC 1521 be- 
gins to approach the reg ime of galaxy groups (> IO^'^Mq: 
Humphrey ct al. 2006| ), but, since it represents one of 
the most luminous galaxies in our sample, this will likely 
represent the upper mass envelope of the EllXr galax- 
ies. While the mass constraints in NGC 1521 are com- 



petitive with the best measurements in Humphrey et al 
( ^006|), we note that t hey are poorer than for NGC 720 
(H 



umphrey et al. |2011), especially for the concentration. 




100 



200 
Radius (kpc) 



500 



Fig. 10. — Baryon fraction profile inferred from our best-fitting 
models (black line). The shaded grey region indicates the l-cr sta- 
tistical uncertainty in the fits. The red shaded region indicates 
the gas fraction profile and l-cr uncertainty. The dashed line in- 
dicates the bpst-fitting Pnsmnlngical va]up of fb , based On the 5- 
year WMAP data ( [Punkley et ah 2009 ) . We indicate the physical 
scales corresponding to Rvir and various other standard radii, and 
the dotted line indicates the radial extent of the data being fitted. 
We note that the fj, profile is relatively flat, and is consistent with 
the Cosmic mean by ^jR,h This is vfiry ymilar to the profiles 



measured in NGC ^ 720 ( Humphrejj | 0t ah . 2011) and the fossil group 
RXJ1159-I-5531 ( Humphrey et al. ' lUTIa T 



tions for NGC 1521 have some role to play in this effect, 
it mostly reflects degeneracies between the dark and lu- 
minous matter, since the dark matter only dominates 
outside ~20 kpc, which is comparable to the scale ra- 
dius (31^i[^ kpc). Nevertheless, the virial mass is more 
tightly constrained, which reflects the good radial cover- 
age at larger scales. 

Even though the light profile of NGC 1521 is quite com- 
plex (§ 2.2), and so a simple interpret ation in terms of the 
toy Se rsic+NFW model proposed by Humphrey fc Buote 
(2010) may not be possible (i.e. the rnodel lines in Fig 
we found that we were able to fit the radial (total) mass 
distribution within lORc with a purely powerlaw den- 
sity profile distribution. In Fig |[ we show the locus of 
NGC 1521 in the a-Rc plane, which lies very close to the 



g-Rp relation established by other systems (Humphrey 



Although the slightly lower fiux and shallower observa- 



& Buote 201C). This further supports the idea that ap 
proximately powerlaw total mass distributions may be a 
natural consequence of the process of galaxy formation. 

6.4. Gas and baryon fraction 

Based on our self-consistent hydrostatic mass model, 
and the fit to the density profile, we were able to con- 
strain the gas and baryon fractions in NGC 1521, allow- 
ing a direct comparison with both NGC 720 and more 
massive groups and clusters. In Fig hoL we show the ra- 
dial profile of fgas and fb, and we tabulate the values at 
various interesting scales. We see qualitatively the same 
behaviour as in NGC 720, i.e. an approximately flat fb 
profile, close to the cosmic mean. When evaluated at 
Rvirj fb is consistent with baryonic closure, in agreement 
with the picture that the majority of the baryons ex- 
ist in the hot ISM, and the behaviour seen in NGC 720 
and also the more massi ve fossil group RXJ1159-I-5531 
(Humphrey et al. 2012a). Within plausible systematic 
uncertainties, it is possible that fb vir in NGC 1521 could 
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Fig. 11. — 1-cr confidence region for the entropy profile 
model of NGC1521, scaled by its characteristic entropy K500 
(=38.5 keV cm^), and shown as a function of Rsoo- We overlay 
deprojected entropy data-points (see § 5.7), similarly scaled. The 
Chandra data are marked with (red) triangles and the XMM data 
with (magenta) stars. The dashed line indicates the "haaeline" prp- 
diction from gravitational structure formation (K^oit et al. 200E). 



Also shown (blue shaded region) is the scaled entropy model, cor 
rected for the gas fraction profile, which agrees better with the 
baseline model. In yollr.i»r alir^.r tlio or,ti-r^py profile for the iso- 
lated galaxy NGC 720 (Humphrey ct al. 2011). 



be slightly lower (~0.10), but even in that case the sys- 
tem has retained the bulk of its baryons. Whether the 
two galaxies are representative of systems at this mass- 
scale, or if they constitute special cases, remains to be 
established. Results for the full EllXr sample will help 
us to address this question. 

Unfortunately, the measured value of fb in NGC 1521 
does not include all of the system's baryons. In par- 
ticular, it ignores the dwarf companions, and a possi- 
ble extended stellar envelope, analogous to intracluster 
light. Without more velocity measurements, it is dif- 
ficult to estimate what fraction of the galaxies within 
the projected virial radius (Fig 0) should actually be in- 
cluded in such a calculation. Nevertheless, for a halo of 
Mvir= 6 X IO^^Mq, anywhere between ~5-40% of the 
total stellar light c ould be in the satel l ite galaxies and 



exten ded envelope ( Purcell et al 
^OOTD , which could 
at most 



2007 



Gonzalez et al 



errors. We note that NGC 152 



ead to an underestimate in fb,vir of 
which is comparable to the statistical 



significant cool gas ( [Huchtmcicr 



itself does not contain 



19941) 



In contrast, fgas is more robustly known, and can easily 
be compared to measurements in other systems. Within 
R-soOj we found that fgas is in good agreement with an 
extrapolation of the trends seen in galaxy groups and 
clust ers: extrap olating the fgas-Rsoo relation from Gio 



dini j;t al.| (|00g) to the mass of NGC 1521, we would ex- 
pect fgas=0.046 ± 0.008, very close to the observed value 
OMOt^^nnt Similar behaviour was seen in NGC 720 



(Humphrey et al. 2011). 



6.5. Entropy profile 

As expected for approximately hydrostatic gas, we find 
that we obtain a good fit to the data with a model re- 
quiring a monotonically rising entropy (S) profile. We 
show the model profile in Fig |ll| (grey shaded region). 



scaled by the "characteristic entropy" i^sooi and shown 
as a function of fraction of R500 reached. We find that 
the characteristic shape of the profile is qualitatively sim- 
ilar to that of NGC 720, showing a central "plateau", 
and gradual steepening of the profile with radius, ap- 
proximately reaching S cx R^'^, although it does not 
exhibit the flattening outside '^O.lRsoo in that galaxy. 
At all radii, the entropy profile is significantly enhanced 
over the "baseline" r nodel for gravity -onlv structure for- 
mation simulations ( Voit et al. |2005 ), indicating signifi- 
cant entropy injection. Given the more massive halo of 
NGC 1521, as compared to NGC 720, it is unsurprising 
that the entropy injection appears to have affected the 
profile more significantly in the latter. 

To provide a less model-dependent view of the entropy 
profile, we overlay in Fig[ll] a series of data-points, which 
are directly computed from the deprojected density and 
temperature profiles. These were obtained by emulating 
the behaviour of the "projct" Xspec model, and correct- 
ing for emission projected into the lin e of sight from out- 
side the outermost annulus (see § ^]^) . These data agree 
well with the smooth model, giving us confidence in our 
projection p rocedure. 

Following [Pratt et al.| ( |201C| ), we investigated whether 
scaling the entropy profile by a correction factor 

{{fg{R)/fb,u)^^^, where fg{R) is the gas fraction pro- 
file, and fb,u is the Cosmic baryon fraction, 0.17) brings 
it into better agreement with the baseline model. We 
show this "fgas-corrected" entropy profile in Fig p], which 
clearly agrees much better with the baseline model. This 
is consistent with a picture in which entropy injection 
primarily manifests itself by pushing the gas out to large 
radii. 

6.6. Abundance profile 

An intriguing result from our study is the detection of a 
negative abundance gradient in the hot ISM of NGC 1521 



grouDS (e.g. 


Humohrcv & Buotc 


200d: Buott b002: 


Buote 


it al. 2003, 


2004; 


Rasmussen & PonmanI 2067(), and the 


~Milkv Wav-mass elliptical NGC 720 (|Humphrey et al. 



be commonplace even in isolated (X-ray bright) ellipti- 
cal galaxies, which has implications for the b ulk gas mo- 



tions responsib le for distributing the metals (Mathews 
Brighenti 2003), and may i mply large-scale flow s driven 
by low-leveTMJN heating ([Mathews et al.[[2004| ) 



In general, the emission-weighted Fe abundance (Zpc 
of the hot ISM in an early-type galaxy is co nsistent with. 



or higher th an, the metallicity of the stars (Humphrey & 
Buotc 2006). In the case of NGC 1521, however, there 
appears some tension between the two measurements; 
using XMM, the Fe abundance, even averaged over the 
central ~10 kpc is only 0.44±0.09 times Solar, while the 



6.2 



results from fitting the published Lick indices (see 
Ogando ct al. 2008) imply a global stellar abundance of 
1.4±0.3 times Solar, making this system an outlier (^3.6- 
<T away from equality between the m etallicitv of the stars 
and ga s) from the relation found by Humphrey fc Buote 
(200™. One possible solution to this problem is the "Fe 
bias" , a systematic underestimate of the Fe abundance 
when a single temperature plasma model is fitted to an 
inherently multi-temperature X-ray spectrum (Buote fc: 
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Fabian 199i; Buote| pOOOb ). Although the temperature 
profile ot iNGC1521 is close to isothermal, we experi- 
mented with adding an additional ~0.2 keV gas compo- 
nent to the central XMM h'm. The fit improved modestly 
(AC=9.7 for 2 degrees of freedom), with ^30% of the gas 
in that bin in the cooler phase, and we found a subtle 
increase in the best-fitting abundance, accompanied by 
a significantly expanded error bar (ZFc=0.50to os)- The 
increase of the error bar sli ghtly eased tension with th e 
other objects discussed in Humphrey &: Buote ( 2006 ), 
making it only ^2.7-a discrepant with parity between 
the metallicities of the gas and stars. We note that, 
although this would suggest non-isothermal gas in the 
central region, one can still obtain reliable mass profiles 
with our modelling procedure by interpreting a single- 
temperature fit to the data, provided care is taken to av- 
erage the model suitabl y, as we did in the present work 



( pastaldello et al.||2007| ). 
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